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In the framework based on the quasipotential method and relativistic quark model 
a new covariant expression for the heavy quark fragmentation amplitude to fragment 
into the pseudoscalar and vector S-wave heavy mesons is obtained. It contains all 
possible relativistic corrections including the terms connected with the transforma- 
tion law of the bound state wave function to the reference frame of the moving 
meson. Relativistic corrections of order p 2 /m 2 to the heavy quark fragmentation 
distributions into (cc), (be) and (bb) states are calculated as functions of the longi- 
tudinal momentum fraction z and the transverse momentum pt relative to the jet 
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lO ■ I. INTRODUCTION 
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The decay and production processes of the bound states with heavy quarks are inves- 
tigated with greater intensity in the last years. The research aims of many experiments 
(ALEPH, DELPHI, SLD, CLEO, Belle, SELEX, LHC-b) are directed on the growth of ex- 
perimental accuracy in the derivation of the static characteristics of heavy hadrons, their 
production and decay rates in different reactions jH Q, 0, 01- The production of heavy 
mesons and baryons via heavy quark fragmentation in the e + e~ - annihilation represents 
one of the possible mechanisms for the formation of heavy hadrons with two heavy quarks. 
The fragmentation cross sections for the production of heavy hadrons can be calculated in an 
analytical form using the factorization hypothesis. The heavy quark production amplitude 
can be calculated on the basis of perturbative QCD. The characteristic quark virtualities 
of heavy quarks in the hard production are of the order of their masses while the quark 
virtualities in the bound state are much less than their masses due to the nonrelativistic 
motion. So, the total amplitude can be represented as a convolution of the hard transition 
amplitude with a nonperturbative factor (the wave function) determining the transition 
of free heavy quarks into a bound state p. The fragmentation mechanism was used for 
the study of the production processes of heavy mesons and baryons in e + e~ annihilation 
in Refs.Jji, 0, HI, 0, EH EH E2 EH (a more complete list of references can be found in 
Refs.0, 0, The growth of theoretical accuracy for the calculation of corresponding 

production cross sections can be reached in two ways. Firstly, it is necessary to take into 
account radiative corrections to the perturbative amplitude describing the production of free 
heavy quarks via heavy quark fragmentation. Secondly, we must consistently consider the 
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relativistic corrections in the fragmentation amplitude connected with the relative motion 
of heavy quarks forming heavy hadron. From the point of view of NRQCD both effects are 
caused by the matrix elements as a function of the typical heavy quark velocity in the bound 
state rest frame of orders 0(vq) and 0(vq), respectively [^]. The experimental data indicate 
that the calculations of different production probabilities for heavy quarkonium and double 
heavy baryons should be improved by a systematic account of relativistic corrections. Such 
effects as the relative motion of heavy quarks forming heavy quarkonia and diquarks, the 
diquark structure effects in the calculation of the fragmentation functions of heavy diquarks 
should be considered. The role of relativistic effects was studied already in the processes of 
c-quark fragmentation into J/^f, r\ c in Ref.|l4| on the basis of the Bethe-Salpeter approach, 
the gluon fragmentation into S-wave quarkonium in Ref . 0] and in the inclusive production 
of polarized J/\& from 6-quark decay in Ref. 0] . The consideration of the intrinsic motion of 
quarks forming the heavy measons can explain the discreapancy between theoretical predic- 



tions and experimental data for the cross sections of the process e + e — > ^r] c [17. 
aim of the present work is to get a systematically improved description of the relativistic 
effects in the processes of the heavy quark fragmentation in the quasipotential approach 
[l9| . Our goal also consists in the calculation of the relativistic corrections in heavy quark 
b and c fragmentation functions into pseudoscalar and vector heavy mesons (Q1Q2) on the 
basis of the relativistic quark model used earlier in the calculation of mass spectra of heavy 
mesons and baryons and their decay rates in different reactions [20L 21 . 22T |. In particular, 



we investigate double distribution probabilities for the heavy quark fragmentation over lon- 
gitudinal meson momentum z and transverse meson momentum p?. Analytical expressions 
for the fragmentation probabilities as functions of transverse momentum of heavy mesons 
(Q1Q2) are obtained. 

II. GENERAL FORMALISM 

The heavy meson production through the heavy quark fragmentation is shown in Fig.l. 
On the first stage Z° boson decays into a quark- ant iquark pair with four-momenta q and q' 
respectively. After that one heavy quark with the four-momentum q fragments to the heavy 
quarkonium. In the quasipotential approach the invariant transition amplitude of a heavy 
quark b or c into a heavy meson can be expressed as a simple convolution of a perturbative 
production amplitude T{y\ , p' 2 , v a ') of free quarks and the quasipotential wave function of 
the bound state (Q1Q2) *p(p) HEB 

M{q,P,p',q') = J ^ 1 |r p (p)T(p' 1 ,p / 2 ,p / ,g / ), (1) 

where four-momenta of fragmenting quarks (6, c) and spectator antiquarks (b, c) forming the 
heavy meson are defined as follows: 

Pi = Vi p + P, P2 = V2P ~ P, (2) 

p' is the four-momentum of a free spectator quark b or c and P is the four-momentum of 
the heavy meson. The coefficients rji^ in the definition (2) are taken in such a way that the 
following orthogonality condition is fulfilled: 

(p • P) = 0, 7& l2 = ' ( 3 ) 
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FIG. 1: The Feynman diagram for the fragmentation of a heavy quark Q\ with a four- momentum 
q to a heavy meson (Q1Q2) with a four-momentum P. T is the vertex function determining the 
production of the quark-antiquark pair in the Z° decay. 



M = (mi + m 2 + W) is the bound state mass. 

The transition of the pair of a heavy quark and antiquark into color-singlet mesons can 
be envisioned as a complicated process in which the colors and spins of the heavy quark 
and antiquark play an important role. Different color-spin nonperturbative factors entering 
the amplitude T(p 1 ,p 2 ,p',q') control the production of the heavy quark bound states. In 
this process the gluon virtuality k 2 ^> Aq CD and the strong coupling constant a s {k 2 ) <C 1. 
Then the hard part of the fragmentation amplitude (1) in the leading order over a s takes 
the form: 



(4) 



3^ (s- m\ 

where Y is the vertex function for the transition of the Z° boson into the quark-antiquark 
pair; the gluon propagator is taken in the axial gauge with four-vector n = (1, 0, 0, —1): 



D\ a {k) = 



k 2 + ie 



k a n x + k\n a 



k ■ n 



(5) 



s = q 2 , k = (q — 7]iP — p) = (i]2P — p + p') is the gluon four momentum. The color factor 
(T a )u(T a ) m jSij I \/3 = 45 m ;/3\/3 was already extracted in the amplitude (4). The transfor- 
mation law of the bound state wave functions from the rest frame to the moving one with 
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four- momenta P is given by 2^, 23| 

^(P) = Dl /2 ' pa (RZ) D 2 2 ' ^(<)<(p), (6) 

*^(p) = * £ o T (P)^ + 1/2 ' £A (^ P )^2 + 1/2 ' TCT « P ), 

where i? w is the Wigner rotation, Lp is the Lorentz boost from the meson rest frame to a 
moving one, and the rotation matrix D 1//2 (i?) is defined by 



'l 
,0 



J)a 1 , / 2 2 «) = 5- 1 ( Pi , 2 )5(P)5( P ) 



(7) 



where 



s(p) = \ e{P l + m fi+ (ap) 



2m y e(p) + m J 

is the usual Lorentz transformation matrix of the four-spinor. For further transformations 
of the amplitude (4) the following relations are useful: 

S aP (A)u*(p)= J2 vr*mD$(RZ), (8) 

(T=±l/2 



U 



(p)S£(A) = E ^aV /2 (<K(Ap) 



<T=±l/2 

Using also the transformation law of the Dirac bispinors to the rest frame 

Mi(p) = Mi(0) ===== Pi = (ei,p), 

v /2ei(p)(e 1 (p)+m 1 ) 

f 2 (- P ) = , ===== v 2 (0), P2 = (e a , -p), (9) 

/2e 2 (p)(e 2 (p) +m 2 ) 

we can introduce the projection operators fl p,v onto the states (Q1Q2) in the meson with 
total spin and 1 as follows: 

TL P > V = [^(0)« 1 (0)]s=o,i = 75(6*)^. (10) 

As a result the heavy quark 6(c) fragmentation amplitude into the mesons (be), (66) or (cc) 
takes the form: 



„w D / a 2a,V2M /■ D A „(fc) dp - , (p 2 - m 2 ) , v 

M(q,P,p,q) = / ? 277^3^° (PJW 7a -7= = x 11 

3V6 V (s-mi (2tt 3 J2e 2 (p)(e(p) +m 2 ) 



xe~>)(6 + 1) . === = 7A(g + mi)r^(g'), 
v /2e 1 (p)(ei(p)+mi) 



where the four-vectors e, pi j2 are given by: 

1 +t»° 



e = L P (0,e)= ev ,e+f^ , (12) 
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P h2 = S(L P )m 2 S- l (L P ), S(L P )(1 ± j^S-^Lp) = ±(v ± 1) 



P 
M' 



Transforming the bispinor contractions in the numerator of the expression (11) we can 
find the following expression for the heavy quark fragmentation amplitude including the 
effects of relative motion of the heavy quarks: 



M(g,P,p',g') 



2a s V2M 



*o(p) 



D X(T (k) dp 



3a/6 J [EM (*i(p)+™i) J EM Mp)+m 2 ) (g _ m j) (2^)3 
V mi 2mi V m2 2m,2 



(13) 



xm 2 (p')7 ( 7 
'v + 1 



v - 1 



+ v 



p 



p 



2m 2 (e 2 

2 



+ 



m 2 J 
V 



(v)(v + l)x 



2m! (ei + mi) 2mi 



2m 2 _ 

jx(q + mi)r Q fi(g'). 



The fragmentation amplitude (13) keeps at least two sources of relativistic corrections. 
The corrections of the first group appear from the quark- ant iquark interaction operator. 
They can be taken into account by means of the numerical solution of the Schroedinger 
equation with the relevant potential. The second part of these corrections is determined by 
several functions depending on the momenta of the relative motion of quarks p. In the limit 
of zero relative momentum p the amplitude M(q,P,p',q') was studied in Ref. BQB. At 
last, there exist the one-loop corrections to the fragmentation amplitude which we have not 
considered here. 

Heavy quarkonium can be characterized by the hard momentum scale m (the mass m of 
the heavy quarks), the soft momentum scale ijivq and the ultrasoft momentum scale mvq. 
We assume that the heavy quarkonium is a nonrelativistic system. This implies that the 
ratio p 2 /m 2 ~ Vq <C 1. So, we introduce the expansion of all factors in Eq.(13) over relative 
momentum p up to terms of the second order: 



1 

T 2 



11 4 



I.' 2 
h 



1 

kn 



1 



+ 



pn 



r] 2 s + r]im\ - r/i?7 2 M 2 , 
(pn) 2 



(14) 



(qn — rjxPn) (qn — r/iPn) 2 (qn — rjiPn) 3 

J2 



(px + mi) 



m\(v + 1) + £77- h p, (p 2 - ra 2 ) = m 2 (v - 1) + v P 



P 



2mi 
L P (0,p) 



2m 5 



P- 



pv,p + 



v(pvj 



V 







Substituting the expansions (14) into Eq.(13) we obtain: 



M(q,P,p',q' 



2a sV / 2M 
3^ 

XU 2 (p')j a 



1 - - p2 



3 p 2 



2 



8 8m 2 , 



dp 



(s - ml) (2tt) 3 



x 



(15) 



1 M* P 

h V 



P 



4mo 



2m 2 



[v)(v+ l)x 



6 



v + 1 



+ 



P 



Am{ 2rrii 



'j x (q + m 1 )T a v 1 (q')x 



x 



k 2 



k 4 1 



L6 
A 



+ 



pn 



+ 



(gn — rjiPn) (qn — rjiPn) 2 (qn — rjiPn) 3 

Let us emphasize that for the system of two heavy quarks relative motion corrections 
entering in the gluon propagator or heavy quark propagators have the same order 0(v%) 
contrary to the system including heavy and light quarks. In the last case leading order 
relativistic corrections are determined by the relativistic factors belonging to the light quark 
but other terms contain an additional small ratio m q /mQ. The obtained relation (15) which 
has the form of a three dimensional integral in the momentum space is valid when the inte- 
gration is restricted to the soft momentum region, where the wave function has a significant 
support. Otherwise it would diverge at high momenta. Moreover, our aim consists in pre- 
serving here only the terms of the second order over |p|/m omitting corrections of higher 
order. 



III. HEAVY QUARK FRAGMENTATION FUNCTIONS INTO P- AND 

V-MESONS 

We use the NRQCD factorization approach to the calculation of the fragmentation reac- 
tions which was developed in Refs.0, 0|. The fragmentation function of the heavy quark Qi 
to produce 1 Sq or 3 Si (QiQ 2 ) meson states is determined by the following expression: 

1 f , „ ( M 2 m\ \ , IMP , . 

where q is the energy of the fragmentating quark: q^ = (q , 0, 0, Jq^ — s) ; M = ■u 1 (g)rt' 1 (g / ) 
is the amplitude of free quark Q 1 production on the mass shell, z is the meson longitudinal 
momentum fraction relative to the fragmenting heavy quark. Let us consider the fragmenta- 
tion production of the vector meson. Omitting the momentum of the relative motion of heavy 
quarks p in Eq.(14) we obtain the fragmentation amplitude which contains the leading order 
contribution and the correction due to the quark bound state energy W (M = mi + rri2 + W): 



= 1 x ( 

3V6 {s - m() (r} 2 s +rj 1 mi- rji^M 2 ) 

u 2 {p')2t{v){q + m l )Y a v 1 {q') + {s + ^^^^ ^V)™^ + 1)^1^) + 

+7 T -rU 2 {p)ne {v){v+l)pT a v 1 {q)+- — -u 2 {p)ne (i;)(u+l)(g+mi)r a i;i(g ) 

(nq — rjxnP) [nq — rjinP) 

Taking into account linear terms in the binding energy W in the expansion of Eq.(16) we 
next perform an averaging and summation over the meson spins in initial and final states 
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in the square modulus \M\ 



1 1 

spin 



(18) 



and then consider the limit go ^ oo in the obtained expression. In this limit P and q have the 
order of the Mz mass and the coefficients in corresponding expressions are determined by the 
heavy quark masses of the fragmenting quark mi and the spectator quark m 2 . In the leading 
order we can substitute P = zq, and the trace in \M\ 2 over the Dirac indices is proportional to 
Tr(T a q'Tpq'). It will disappear in the ratio |M| 2 /|M | 2 . Then the fragmentation probability 
(16) can be written as a sum of two terms: 



D 



Qi-+V(QiQ 2 ) 



8a 2 |^(0)| 2 rz(l-zf 



(19) 



27m| [1 - (1 - r)zf 

v = 2-2(3-2r)^ + 3(3-2r + 4r 2 )2 2 -2(l-r)(4-r + 2r 2 )z 3 + (l-r) 2 (3-2r + 2r 2 )z 4 , (20) 



W 



Vl 



3m2[l — (1 — r)z} 2 



-12+6r+(60-84r+48r 2 )^+(-138+259r-296r 2 +102r 3 )z 2 + (21) 

+ (192 - 420r + 552r 2 - 536r 3 + 330r 4 )^ 3 + (-168 + 446r - 614r 2 + 746r 3 - 612r 4 + 232r 5 )z 4 + 

+ (84 - 274r + 476r 2 - 650r 3 + 574r 4 - 272r 5 + 62r 6 )z 5 + (-18 + 57r - 126r 2 + 282r 3 - 

-412r 4 + 333r 5 - 144r 6 + 28r 7 )z 6 + r(10 - 40r + 56r 2 - 24r 3 - 14r 4 + 16r 5 - Ar 6 )z 7 . 

As mentioned above there exist several sources of relativistic corrections |p| 2 / m i,2 m the 
expression (15). The first part of terms appears from the expansion of the gluon propagator 
and relativistic factors in the Dirac bispinors. The structure of the spinor matrix element 
in this case is the same as in Eq.(17): 



2a s \[2M 



dp (mi + m 2 ) 



x 



M2l(g ' P ' pW) = "W J ^^ M (^rn 2 (s-rn 2 ) 2 
A(qpf p 2 p 



X 



(22) 



8m 2 8m 2 



+ 



u 2 (p')2~e (v)(q + m 1 )T a v l (q') + 
2(qp)(pn) (pn) 2 



! lf_ 4(gp) 2 _ jr_ _ jr_ _ _ 

kg kf> § m i § m 2 k 2 (nq — rjinP) (nq — n\nP) 2 



x 



x ■ 



v s — m\) 



-u 2 (p')ne (v)(v + l)I>i(g')}- 



(nq — rjinP) 

Another part of corrections is determined both by the gluon propagator terms and relativistic 
addenda p in the square brackets of Eq.(15): 



M 22 (q,P,p',q') = 



x 



(np) 



+ 



2qp 



(nq — rjinP) k% 



2a s y/2M 

(Panx + p x n a 
(nq — rjinP) 



- dp (mi + m 2 ) 

^ (0lQ2) ' o(P) (2vr)3m 2 ( S -m 2 ) 2 X 



■u 2 (p')^*(v)(v + l)7 A (g + mi)r a ui(g'), 



(23) 
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2a s \/2M 



M 23 (q } P,p', q ') = |^=- J * (0lQ2 ),o(P 



(PaTix +P\n a ) _ , ,* 

x -7 w ha 

(rig — r\\nr) 

M 24 (q,P,p',q')-- 



'' T{v){v + l) + e{v){v + l) '' 



2m 2 

2a s V2M 
3^/6 



dp (mi + m 2 ) 
(27r) 3 m 2 (s — mf ) 2 

7a (g + mi)r a 7Ji(g') 



(24) 



2mi 



- <ip (mi + m 2 ) 



(27r) 3 m 2 (s — m 2 ) 2 



(25) 



x 



9\a 



(k a n x + k x n a ) 
(nq — rjinP) 



U2{p')lo^-e* {v){v + x )^7A(g + mi)r a 7Ji(g')- 



Further transformations of the interference terms [M 2 \M 2i + Mfj M*^ ) (i = 2,3,4) in the 



square modulus \M\ are performed by means of the system Reduce The relation 



P 2 - 3P° 2 

9P° 2 



(26) 



which accounts the orthogonality condition (3) is used for the integration over the angle 
variables in the relative momentum space. Omitting numerous intermediate expressions 
appearing in the calculation |M| 2 which have sufficiently cumbersome forms we present the 
final result for the relativistic correction in the fragmentation function (16): 



j~\rel 



iQ 2 )( 2: ) 



8a 2 |*(0)| 2 rz(l-zf 



• v 2 , 



(27) 



v 2 (z) 



(P 2 ) 



27ml [l-(l-r)z] 6 ' 
-2 + 16r (2 - 3r) + 2(1 - 2r) (5 - 44r + 30r 2 )z + (28) 



36[1- (l-r)^ 2 (l-r) 2 m| 

+ (-23 + 18r - 206r 2 + 508r 3 - 432r 4 )z 2 + 8(1 - r)(4 + 40r - llr 2 - 50r 3 + 66r 4 )^ 3 - 

-2(l-r) 2 (14 + 214r + 47r 2 -248r 3 + 120r 4 )2 4 + 2(l-r) 3 (7 + 99r+120r 2 -230r 3 + 108r 4 )^ 5 - 

-(1 - r) 4 (3 + 26r + 82r 2 - 84r 3 + 48r 4 )z 6 ] , 

here (p 2 ) is the special parameter determining the numerical value of relativistic effects (see 
the discussion before Eq.(60)). Integrating expressions (19), (27) over z we obtain the total 
fragmentation probability: 



- D Qi->V(QiQ 2 )( Z ) dz 
•J u 



8ag|tt(0)| 2 
405m 3 ,(l - r) 6 



' Mr) + —fiv(r) + { -^f 2v (r)) , (29) 
m 2 mi I 



f 0v (r) = 24 + 85r - 235r 2 + 300r 3 - 85r 4 - 89r 5 + 15r(7 - 4r + 3r 2 + 10r 3 + 2r 4 ) In r, (30) 
Mr) 1 



42(1 - r 



2607 + 7185r + 23576r 2 - 1 16018r 3 + 159670r 4 - 170373r 5 + 153860r 6 - 



(31) 



-34906r 7 - 20387r 8 + 210r(-42 + 223r - 388r 2 + 236r 3 - 268r 4 - 43r 5 + 212r 6 + 28r 7 ) lnr 



Mr) 



1 



252(1 -r) 2 



488-3017r+48979r 2 -201740r 3 +136955r 4 +23597r 5 -20958r 6 +15696r 7 - 

(32) 
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-105r(7 - 200r + 347r 2 + 1194r 3 - 222r 4 + 156r 5 + 48r 6 ) lnr 

The parameter r = 777,2/ (7711 + 7772) is the ratio of the constituent mass of a spectator quark to 
the mass of two quarks composing the heavy meson. The contributions to the fragmentation 
functions (20)-(21) and (28) give the distributions in the longitudinal momentum z of the 
heavy mesons. They are shown in Fig. 2. We find it convenient to present here the fragmen- 
tation functions of antiquarks for the comparison with experimental data. To extend the 
present calculations to the distributions in the transverse momentum px of the the heavy 
meson, the following relation between the invariant mass s of the fragmenting heavy quark, 
the transverse momentum pt and longitudinal momentum z is taken into account 25] : 



z 1 — z 



(33) 



Introducing further the dimensionless variable t = pr/( m i + m 2)i we can determine the 
Pr-dependent relativistic corrections to the fragmentation probability •Dq 1 -^v(q 1 q 2 ) (t) j^ : 



iQ 2 )(*) 



4a 2 I 



dz- 



2MH 
zT^z 



D Qi^V(Q 1 Q 2 ) {Z,s(z,t)) 



(34) 



27m 3 ,(l -r)H 



W (p 2 ) 
D 0v (t) + —D lv (t) + ^-D 2v (t) 

7772 7772 



D 0v (t) = 
+ arctan 



-30r 3 +30r 4 )/j+(61r-33r 2 -48r 3 +20r 4 )t 3 + (-5+13r-16r 2 +4r 3 +4r 4 )/; 5 + (35) 
t(l-r) 



30r 4 + (-99r 2 - 66r 3 + 30r 4 )t 2 + (9 + 20r + r 2 + 22r 3 + 8r 4 )t 4 + 

+ 



r + t 2 

+ (9 - 12r + 4r 2 + 8r 3 )t 6 l + ln(r) [-96r 3 t + (48r + 56r 2 + 16r 3 )t 3 



+ ln- 



;i + t 2 



r 2 + t 2) 



-48rH + (24r + 28r 2 + 8r 3 )t 3 + (-4 + 8r - 4r 2 )t 7 



D lv {t) 



r 24rt( 1 — r) 6 

{ ^\ - ' +(24r-144r 2 + 360r 3 -480r 4 + 45r 5 + 486r 6 -291r 7 )t + 



3i 2( 1 _ r )2L 1 + f 2 

+ (-24r + 144r 2 - 633r 3 + 2115r 4 - 1233r 5 - 1266r 6 + 897r 7 )t 3 + 
+ (-18r + 1130r 2 - 4129r 3 + 3938r 4 + 595r 5 - 1994r 6 + 478r 7 )t 5 + 
+ (30 - 215r + 581r 2 - 742r 3 + 738r 4 - 596r 5 + 148r 6 + 56r 7 )t 7 + 



(36) 



(60r - 255r 2 + 360r 3 - 195r 4 + 30r 5 )rj 9 



+ 



+ arctan • 



t(l-r) 



, 315r 6 - 315r 7 + (90r 4 - 1650r 5 - 345r 6 + 960r 7 )t 2 + 
(r + t z ) L 

+ (-63r 2 - 1800r 3 + 5646r 4 - 714r 5 - 3495r 6 + 741r 7 )t 4 + 
+ (-54 + 237r - 33r 2 - 807r 3 + 642r 4 - 621r 5 + 423r 6 + 150r 7 )t 6 + 
+ (-54 + 297r - 507r 2 + 216r 3 - 252r 4 + 147r 5 + 126r 6 )t 8 + 



4\ .10 



+ (-60r + 195r 2 - 165r 3 + 30r 4 )t 



+ 
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+ ln(r) 



'-576r 5 + 288r 7 )t + (-432r 3 + 48r 4 + 5856r 5 - 4416r 6 - 48r 7 )t 3 + 



+(-144r + 1296r 2 - 2136r 3 - 768r 4 + 1320r 5 + 432r 6 )t 5 



+ 



+ ln 



r 2 + t 2 ) 



-288r 5 + U4r 7 )t + (-216r 3 + 24r 4 + 2928r 5 - 2208r 6 - 24r 7 )t 3 + 

+ (-72r + 648r 2 - 1068r 3 - 384r 4 + 660r 5 + 216r 6 )t 5 + 
+ (24 - 156r + 312r 2 - 192r 3 + 72r 4 - 60r 5 )t 9 ] }, 

D 2v (t) = r— -|-420r 5 (l - r)t + (6074r 3 - 20646r 4 + 11976r 5 + 2596r 6 )t 3 + (37) 

36^(1 — r) 2 1 

+ (-573r + 2309r 2 + 5008r 3 - 5180r 4 - 1636r 5 + 72r 6 )t 5 + 
+ (5 - 189r + 1124r 2 - 960r 3 + 68r 4 + 48r 5 - 96r 6 )t 7 + 
t(l -r) 



+ arctan ^ 420r 6 + (-7470r 4 + 15480r 5 + 5040r 6 )t 2 + 

r + 1 2 1 

+ (1887r 2 - 1758r 3 - 20334r 4 - 6936r 5 + 900r 6 )t 4 + 
+ (-9 + lOOr + 1303r 2 - 1614r 3 - 632r 4 + 216r 5 - 168r 6 )t 6 + 



-9 + 108r - 1336r 2 - 96r 3 + 64r 4 - 192r 5 )t 8 



+ 



+ ln(r) (-3072r 5 + 4608r b )t + (6240r 3 - 10688r 4 - 24256r 5 + 768r b )r+ 



+ -240r + 712r 2 + 5232r d + 3392r 4 + 384r° )t 



5\ + 5 



+ 



+ ln 



(r 2 + t 2 ) 
{t 2 + 1) 



-1536r 5 + 2304r b )t + (3120r 3 - 5344r 4 - 12128r 5 + 384r b )r+ 



+(-120r + 356r 2 + 2616r 3 + 1696r 4 + 192r 5 )t 5 + 

+ (4 + 72r + 324r 2 - 208r 3 + 96r 4 )£ 9 ] }. 

The leading order contribution (35) coincides with the result of Ref.|26j| and the two other 
terms (36) and (37) determine the relativistic and bound state corrections. The functions 
(35)-(37) are plotted in Fig. 4. 

Let us next consider the calculation of the fragmentation functions into pseudoscalar 
heavy mesons rj c , B c , rjb. The general expression of the fragmentation amplitude consists of 
three terms: 



2V2Ma,|^(0)| (m 1 + m 2 ) 



3VQ m 2 {s-m 2 1 ) 2 
+a 2 u 2 {p')^T a vi{q') + a 3 u 2 (p')j 5 n(q + mx)Y a v 1 {q') 



(38) 



(the fourth possible term u 2 (p')-y 5 nT a Vi(q') does not contribute to the fragmentation func- 
tion) where the coefficients Oj (i=l,2,3) contain the leading order contribution and corrections 
proportional to (p 2 ) and W: 



ai = 1 — 



(P 2 ) 



+ 



(P 2 )r5 , 2 



9 (s — m 2 )r m 2 



r r + 

24 9 9 



(39) 
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1 o 2 2 1 2 1 7 2 2 

l-s(l-r) V2^ "6 r -18 r " + 18 r + r-9 



+ 



2(mi +m 2 )W^ W / 1 z 2rz 

-(l-r) + - - 2 



(s — mf ) (mi + m 2 ) y r 1 — z(l — r) 1 — z(l — r) 

2\ / 1 c 17\ O /^2\ ^2^2 



P 2 ) f 1 2 , 5 17\ , 2(p 2 ) 



r z 



mf V 6 18 12/ 9 [1 — z(l— r)J J 

(P 2 ) 1 / 2 2 2 2 1 1 1 IN 

H n : r — r z H — r H — r z -\ — + 

m? l-z(l-r) V 9 9 3 3 9 9/ 

, (P 2 ) rn 2 2 _ J_ + 1 + J L + 1 _ _ 4 <P 2 > , 



m\ (s — ml) \9 9rz r 9r 2 ^ 9r 2 3^ 3/ 9r (s — m 2 ) 

+ 2(,n 1+ m2 )W-- W_/ 2 _1 

(s — mf ) (mi + m 2 ) \ r 

^ = 1 Wfi\ + W__^_^.W_ (41) 

9mf \ 8/ mim 2 1 - 2(1 - r) 9 
2 (p 2 ) | 2(m 1+ m 2 )^ (1 | ^ / 1 



9r (s — mf) (s — m 2 ) (mi + m 2 ) V r, 

The exact expressions (39)- (41) are obtained on the basis of Eq.(15) (e — > 75) keeping the 
terms of the second order over relative momentum p and binding energy corrections. The 
squared modulus amplitude IM3I 2 leads by using the definition (16) after the integration 
over s to the following fragmentation distribution for the production of the pseudoscalar 
heavy mesons: 

, , 8a 2 |^(0)| 2 rz(\ - zf . 

p = 6 + 18(2r-l)^ + (21-74r + 68r 2 )z 2 -2(l-r)(6-19r + 18r 2 )z 3 + 3(l-r) 2 (l-2r + 2r 2 )z 4 , 

(43) 

Pl = (_i2 + 6r + 60z + 24r(-5 + 2r)2+[-126 + r(425-482r + 218r 2 )],2 2 + 

m 2 [l — (1 — r)z\ 2 1 

(44) 

+2[72+r(-329+2r(298+r(-283+131r)))]2 3 -2(l-r)[48+r(-219+r(410+7r(-61+33r))]^ 4 
+2(l-r) 2 [18+r(-79+2r(69+13r(-5+3r)))]z 5 +(r-l) 3 [6+r(-25+r(35+6r(-4+3r)))]z 6 }, 

k = aed-r^-d-rM' i- 6 - 48r - 6| - 5 + 2r{ - 3 + r{ - 3 + 22r))|: - (45) 

-[63+2r(-191+r(431+4r(-197+146r)))]z 2 +8(l-r)[9+r(-108+r(217+5r(-42+29r)))]2 3 - 

-2(1 - r) 2 [24 + r(-390 + r(697 + 4r(-53 + 5r)))]z 4 + 2(r - l) 3 [-9 + r(155 + 2rx 

(-122 + r(-7 + 36r)))]z 5 + (1 - r) 4 (-3 + 42r - 58r 2 + 24r 4 )/}. 

Integrating expressions (43)- (45) over z we obtain the total fragmentation probability for 
the pseudoscalar mesons: 

n P = I iJ^wM* = 405 4 ( ^' r)a (Mr) + -Mr) + ^Mr>) , (46) 
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/pp(r) = 8 + 5r + 215r 2 - 440r 3 + 265r 4 - 53r 5 + 15r(l + 8r + r 2 - 6r 3 + 2r 4 ) lnr, (47) 

/ lp (r) = — [897 - 3640r + 8400r 2 - 59850r 3 + 147105r 4 - 132762r 5 + 46830r 6 - (48) 

42(1 — r) 

-6980r 7 + 210r(6 - 3r - 129r 2 + 115r 3 + 123r 4 - 84r 5 + 18r 6 ) lnr 

U Jr) = -77 [614 - 9051r - 9681r 2 + 106470r 3 - 73815r 4 - 45129r 5 + (49) 

756(1 — r) 2 L 

+41146r 6 - 10554r 7 + 105r(-3 - 204r + 281r 2 + 798r 3 - 402r 4 + 48r 5 + 24r 6 ) lnr 

The transverse momentum fragmentation functions for the production of pseudoscalar 
mesons r] c , B c and r)b can be derived in a similar way as for the vector mesons. The corre- 
sponding expressions are given as follows: 

D Qi^P(.QiQ2)(t) = J dz z ^ _ z ^ D Qi^P(QiQ2) ( z , *)) = ( 50 ) 



4a 2 |#(0)| 2 r 



W (p 2 ) 
D 0p (t) + —D lp (t) + ^i-D 2p (t) 
m 2 wij 



81m 3 (l -r) 6 t 6 

Dpp(t) = {-(1 - r)t[30r 3 - r(61 - 20r + 28r 2 )t 2 - (3 - 48r + 48r 2 - 12r 3 )t 4 ] + (51) 

^ ~ T ) [lOr 4 - 3r 2 (ll + 2r + 2r 2 )t 2 + (3 + 4r + 19r 2 - 6r 3 )t 4 + 
r + t 2 L 

+ (3 + 12r - 20r 2 + 8r 3 )t 6 )l - 24rt ln(r) [4r 2 - (2 + r + 2r 2 )t 2 ]- 



+3 arctan • 



-12tln (1 2 +t2 2 } [4r 3 - r(2 + r + 2r 2 )t 2 + (1 - r)H 6 ]}, 



(r 2 + t 2 ) 

D lp (t) = 2/ / { 24rt(1 ; r)4(2r ~ 1) +3r(8-48r + 112r 2 -128r 3 + 37r 4 + 19r 5 )t+ (52) 
t z (Y — r) 1 1 + t z 

+r(-24 + 144r - 691r 2 + 1032r 3 + 84r 4 - 545r 5 )t 3 + 

+r(10 + 360r - 1180r 2 + 1037r 3 - 495r 4 + 268r 5 )t 5 - 

-3(2 - 19r + 117r 2 - 300r 3 + 320r 4 - 140r 5 + 20r 6 )t 7 + 

35r 6 - 10r 4 (-13 + 15r + 25r 2 )t 2 + 



t(l-r) 
+3 arctan ■ 



(r + t 2 ) 

+r 2 (-15 - 281r + 501r 2 + 96r 3 + 79r 4 )t 4 + 
+ (-6 + 31r - 7r 2 + 79r 3 - 245r 4 + 42r 5 + 8r 6 )t 6 - 
-(6 - llr - 68r 2 + 208r 3 - 160r 4 + 40r 5 )t 8 + 
+24rt ln(r)[-16r 5 + 2r 2 (-7 - lOr + 38r 2 + 4r 3 )t 2 + (-2 + 20r - 21r 2 - 5r 3 - 22r 4 )t 4 ] + 

(r 2 + t 2 ) 

+12* In \ / [16r 6 - 2r 3 (-7 - lOr + 38r 2 + 4r 3 )* 2 + 
1 + t 2 

+r(2 - 20r + 21r 2 + 5r 3 + 22r 4 )* 4 - (1 - r) 2 (2 - 5r + 5r 2 )* 8 ]}, 
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D 2p {t) 



1 



: {-420r 5 (l - r)t + 2r 3 (1957 - 7323r + 4368r 2 + 998r 3 )t 3 



36t 2 (l -r) 2 

-r(429 - 2629r + 152r 2 - 396r 3 + 1684r 4 + 760r 5 )t 5 + 



(53) 



+3(-l + 65r - 308r 2 + 360r 3 



+ arctan 



t(l-r) 



(420r e 



108r 4 - 24r 5 + 16r 6 )t 7 + 



'-5310r 4 + 11640r° + 4440r 6 )t 2 + 



r + t 2 

+ (1023r 2 - 4158r 3 - 9486r 4 - 4080r 5 - 1764r 6 )t 4 + 
+ (-9 - 324r - 9r 2 + 1674r 3 + 1632r 4 - 360r 5 + 72r 6 )t 6 - 
+ (-9 - 36r + 744r 2 - 552r 3 - 24r 4 + 96r 5 )t 8 

+24rtln(r 



-128r 4 + 192r 5 



4r 2 (-35 + HOr + 110r 2 + 26r 3 )£ 2 + 



+ (-10 - r + 150r 2 + 152r 3 + 100r 4 + 4r 5 )t 4 



+ 



+12* In 



(1 + ^ 2 ) 
(r 2 + t 2 ) 



-128r 5 + 192r s - 4r 3 (-35 + HOr + 110r 2 + 26r 3 )* 2 + 



+r(-10 - r + 150r 2 + 152r 3 + 100r 4 + 4r 5 )* 4 - (-1 + Ur - 13r 2 - 4r 3 + r 4 )* 8 )] }. 



The leading order function (51) coincides with the result of Ref. [26j, and the two other 
functions (52) and (53) determine the relativistic and binding energy corrections. Functions 
(51)-(53) are plotted in Fig.5. The obtained results for the transverse momentum distribu- 
tions are valid for the transverse momentum p? up to values of order of the meson mass. The 
further growth of the momentum pt demands the consideration of the omitted corrections 
of order O (pr / Mz ) . 



IV. QUASIPOTENTIAL QUARK MODEL 

To estimate numerical values of the investigated effects in the heavy quark fragmentation 
we used the relativistic quark model. In the quasipotential approach the bound state of a 
quark and antiquark is described by the Schrodinger type equation J2?| 

fb 2 (M) p 2 \ r dq 



\ 2fi R 2/j, R J J (2tt) 

where the relativistic reduced mass is 



*o(p) = / 7 I F(p,q,M)^ o (q), (54) 



E X E 2 M 4 -(m?-w|) 2 

^ " e^Te~ 2 ~ w ' (55) 



and the particle energies Ei, E 2 are given by 



M 2 — m 2 + ml M 2 — m\ + m\ 

1= 2M ' 2= 2M ' ^ ' 
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FIG. 2: The contributions to the fragmentation functions D(c — > J/ty)(z), D(6 — > B*)(z), D(h — > 
T)(z). The thick solid line shows the total fragmentation function, the dashed line shows the 
relativistic correction (27), the dotdashed line shows the correction proportional to the binding 
energy W (21). The thin solid line corresponds to the distributions without corrections. All 
functions have been multiplied by a factor 10 4 . 
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FIG. 3: The contributions to the fragmentation functions D(c — > r) c )(z), D(b — > B c )(z), D(b — ► 
rjb)(z). The thick solid line shows the total fragmentation function, the dashed line shows the 
relativistic correction (45), the dotdashed line shows the correction proportional to the binding 
energy W (44). The thin solid line corresponds to the distributions without corrections. All 
functions have been multiplied by a factor 10 4 . 
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TABLE I: Basic parameters of the relativistic quark model and the total fragmentation probabil- 
ities Q, for the reactions c — ► rj c , J/^f, b — > B c , B*,b^ rjf,, T. 



n 2S+1 Lj 


Particle 


Mass GeF 
[37] 




Round statp 
energy, GeV 


V(0) GeV 3 / 2 


(r> 2 ) GeV 2 


Q 


1% 




2.980 


0.451 


-0.120 


0.27 


0.5 


1.49 x 10" 4 


l 3 ,Si 




3.097 


0.451 


-0.003 


0.27 


0.5 


1.28 x 10" 4 




5 C 


6.270 a 


0.361 


-0.160 


0.33 


0.7 


3.49 x 10" 4 


l 3 5i 




6.332 a 


0.361 


-0.098 


0.33 


0.7 


4.89 x 10~ 4 




Vb 


9.400 a 


0.267 


-0.360 


0.46 


1.4 


0.14 x 10" 4 




T 


9.460 


0.267 


-0.300 


0.46 


1.4 


0.15 x 10~ 4 



"This value was obtained in the Ref. 



here M = E\ + E2 is the bound state mass, are the masses of heavy quarks (Qi and 
Q2) which form the meson, and p is their relative momentum. In the center of mass system 
the relative momentum squared on mass shell reads 



b 2 (M) 



[M 2 - (mi + m 2 ) 2 ] [M 2 - (m 1 - m 2 ) 2 ] 
AM 2 



(57) 



The kernel V(p, q, M) in Eq. (|54|) is the quasipotential operator of the quark- ant iquark 
interaction. Within an effective field theory (NRQCD) the quark-antiquark potential was 
constructed in Ref. jH,|2^| by the perturbation theory improved by the renormalization group 
resummation of large logarithms. In the quasipotential quark model the kernal V(p, q, M) 
is constructed phenomenologically with the help of the off-mass-shell scattering amplitude, 
projected onto the positive energy states. The heavy quark-antiquark potential with the 
account of retardation effects and the one loop radiative corrections can be presented in the 
form of a sum of sp in-independent and spin-dependent parts. Explicit expressions for it are 



given in Refs. |2lL 22| . Taking into account the accuracy of the calculation of relativistic 
corrections to the fragmentation probabilities, we can use for the description of the bound 
system (Q1Q2) the following simplified interaction operator in the coordinate representation: 



V(r) 



4av(jU 2 ) 
3 r 



+ Ar + B, 



where the parameters of the linear potential A = 0.18 GeV 2 , B = —0.16 GeV, 

tV(/i 2 ) = ot s {n 2 ) 



I 4 2 ) n 



(58) 



(59) 



10 



Oi 



T ~ "7T n /> A) = 11 ~ o n f- 



31 

T 9 J ' r " ~~ 3 
Here rif = 3 is the number of flavors and \x is a renormalization scale. All the parameters of 
the model like quark masses, parameters of the linear confining potential A and B, mixing 
coefficient e and anomalous chromomagnetic quark moment k entering in the quasipoten- 
tial V(p, q, M) were fixed from the analysis of heavy quarkonium masses ji^, 21, %| and 
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radiative decays |21[. The heavy quark masses m b = 4.88 GeV, m c = 1.55 GeV and the 
parameters of the linear potential A = 0.18 GeV 2 and B = —0.16 GeV have standard values 
of the quark models. Solving the Schrodinger-like quasipotential equation we obtain an ini- 
tial expression for the bound state wave functions in the case of (cc) , (be) and (bb) systems. 
For numerical estimations of relativistic effects in the production of heavy mesons via heavy 
quark fragmentation we need the values of the wave functions at the origin, the bound state 
energy and the parameter of relativistic effects: / p 2l ^o(p)^p/(27r) 3 . Note that this integral 
would diverge in the high momentum region. Obviously, the reason of this divergence is 
connected with the used expansion of the integral function in the basic equation (13) over 
the ratio p 2 /m 2 . In the coordinate representation this divergence would appear when we set 
r = in the Coulomb part of the potential (58). Different regularization prescriptions are 
commonly used in this case 3(1 3l|, |32[ EH . An approach to the calculation of this integral 



was formulated in Ref. |24| for solving the orthopositronium decay problem in quantum elec- 
trodynamics. Their prescription is in an agreement with the calculations carried out in the 
effective field theories |35|. Unfortunately, in the investigation of the bound states of heavy 
quarks we cannot use it because the valid wave function asymptotics at p — > oo is not the 
Coulomb-like. So, to fix the value of the relativistic correction we explore the dimensional 
regularization scheme where the scaleless momentum integral JV(p — q)^(q) -^Ss re- 



lated to the problem vanishes [30, 31, 36]. Then we can express the necessary quantity in 
the form: 

( p2 > = W) I (^ p2 ^ o(p) = 2mW + 2mlBl m 

The solutions of the Schrodinger-like equation (54) with the potential (58) determine the 
energy spectrum W of the heavy quark system and lead to the numerical values of the 
parameter (60) for the bound states (cc), (bb) and (be) which are presented in Table I. They 
are in qualitative agreement with the other possible approach for the estimation of the value 
(60) based on the natural regularization directly connected with the relativistic structure 
factors entering in the fragmentation amplitude (13). Heavy quark symmetry predicts that 
the wave functions of the vector and pseudoscalar states are different due to corrections of 
order Vq. The analogous statement is valid for the parameter (p 2 ). Nevertheless, in this 
study we neglect this difference and write in Table I equal values for \l/(0) and (p 2 ) for V- 
and P-mesons. Our value (p 2 ) = 0.5 GeV 2 for (cc)-states is slightly smaller than (p 2 ) = 
0.7 GeV 2 used in Ref. El where it was fixed from the analysis of the quarkonium decay 



rates. The theoretical uncertainty of the obtained values (p 2 ) in the Table I is determined 



by perturbative and nonperturbative corrections to the quasipotential [20, 21] and is not 
exceeding 30%. 



V. DISCUSSION AND CONCLUSIONS 

As mentioned above, the problem of heavy hadron production in e + e~ and pp collisions 
became very urgent in last years. The experimental investigations carried out in this field 
allowed to measure the b quark fragmentation function in Z° decays 3]. The study of the 
fragmentation processes is important as a tool to reveal the features of nonperturbative 
quantum chromodynamics. There appear experimental data indicating essential differences 
between the theoretical predictions and experiment 0, 0, 0, E3, 3- In the present study 



we investigated the role of relativistic and bound state corrections in the heavy quark b, c 
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FIG. 4: The contributions to the transverse distributions D(c — > J D(b — ► B*)(t), 
D(b — > T)(t) relative to the heavy quark fragmentation axis. The thick solid line shows the 
total fragmentation function, the dotdashed line shows relativistic correction (37), the dashed line 
shows the correction proportional to binding energy W (36). The thin solid line corresponds to 
the distributions without corrections. All functions have been multiplied by a factor 10 4 . 
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FIG. 5: The contributions to the transverse distributions D(c — > r/ c )(t), D(b — ► B c )(t), D(b — > 
relative to the heavy quark fragmentation axis. The thick solid line shows the total fragmentation 
function, the dotdashed line shows relativistic correction (53), the dashed line shows the correction 
proportional to binding energy W (52). The thin solid line corresponds to the distributions without 
corrections. All functions have been multiplied by a factor 10 4 . 
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FIG. 6: The total fragmentation functions for the production of vector and pseudoscalar heavy 
mesons as a function of z for fi = mi + 2m,2 (dashed lines) and fi = Mz/2 (solid lines). The thick 
and thin solid lines correspond to the vector and pseudoscalar mesons respectively. All functions 
have been multiplied by a factor 10 4 . 
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fragmentation processes. The amplitude of heavy quark fragmentation is obtained in a new 
form (13) which accounts all possible relativistic factors for the calculation of relativistic 
corrections to the fragmentation functions. Let us summarize several peculiarities related 
to the calculation performed above. 

1. We obtain the heavy quark fragmentation functions for both heavy quarks b and c 
which fragment to pseudoscalar and vector heavy mesons starting with the meson production 
amplitude (1). 

2. All possible sources of relativistic corrections including the transformation factors for 
the two quark bound state wave function have been taken into account. 

3. We investigated the role of relativistic effects in the fragmentation probabilities over 
two variables: the longitudinal momentum z and transverse momentum px of the heavy 
meson. 

Analyzing the obtained analytical expressions for the fragmentation functions both in 
longitudinal and transverse momentum we can point out that the calculated corrections for 
all vector and pseudoscalar mesons are not exceeding 20 % the leading order contribution. 
The numerical value of the binding energy correction is dependent on the initial choice of 
the heavy quark masses because W = (M — ni\ — 777.2). So, for example, the binding energy 
corrections are extremely small for the charmonium production in our model where Wj/y = 
—0.003 Gev. The relativistic correction is essentially more important in this case. Earlier 
the binding energy and relativistic corrections were studied in Ref.Q for the fragmentation 
functions of a charm quark to decay into T] c and J/^f. The comparison of our results with 
the calculation in Ref . jl4j shows that in the case of the reaction c — > J/ty the binding 
corrections are numerically close, but relativistic corrections connected with the expression 
(60) are essentially different both in the sign and numerical value. Our relativistic correction 
coincides in the sign with the leading order contribution and is numerically three times 
smaller than in Ref.jl^]. Moreover, our numerical estimations of analytical relations are 
based on a different numerical value for the expression (60). In the paper [lj] the parameter 
(p 2 ) was fixed (using the mass of c-quark m c = 1.43 GeV) by means of a condition analogous 
to our equation (60) without the addendum proportional to the parameter B entering in 
the confinement part of the potential. The results obtained in the present study evidently 
show that the relativistic plus bound state corrections lead to the systematic increase of the 
fragmentation probabilities for the pseudoscalar and vector mesons. Our total fragmentation 
functions for the decays c — > (cc), b — > (be), b — > (bb) retain the initial shape of the leading 
order contribution. In the production of vector mesons both corrections proportional to W 
and (p 2 ) have the same sign giving us a more essential modification of the leading order 
contribution in the comparison with the pseudoscalar meson production. 

The fragmentation functions (19), (27), (42) depend not only on z but also on the fac- 
torization scale fi. They should be considered at a scale \x of the order of the heavy quark 
masses. The evolution of the fragmentation functions to the scale \x = Mz/2 is determined 
by the DGLAP equation @: 

^ 2 j^D Q ^ H (z^ 2 ) = jf ^-Pq^q fe^l D Q^H(y,V 2 ), (61) 



where Pq^q(x) is the quark splitting function [39(. The modification of the z-shape of the 
fragmentation functions is shown in Fig.6. The average values of the momentum fraction 
for the production of different heavy mesons at the scale \x = Mz/2 are the following: 
< z>= 0.50 (J/tf), <z>= 0.46 (rj e ), < z >= 0.63 (B*), < z >= 0.59 (B c ), < z >= 0.56 
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(T), < z >= 0.52 ( Vb ). 

In the case of (cc) or (bb) mesons Eqs.(29) and (46) acquire a more simple form: 



n 



v 







W 

+ — ( 134 In 2 

m 2 



V(QQ)( z )dz 

78149 



840 



+ 



32a 2 |^(0)| 2 rll89 
27ml ^ 30 

(p 2 ) /1078 



571n2 + 



(62) 



n, 



D f 



*P(QQ) 



(z)dz 



2 
2 

8a 2 1 * 



78416 \ 

In 2 

945 J 



2\ 



1546 



444 In 2 + 



(63) 



W ( , 15581N p 2 /16, 139\ 
+ — -104 In 2 + +^^-7^ —In 2 

m 2 V 240 / ml \ 9 48 / 



Numerical values of the total fragmentation probabilities are presented in Table I. The evolu- 
tion conserves the integral probabilities tty and Qp of the fragmentation. Using expressions 
(29) and (46) we obtain the significant experimental ratio 



r](r) 



(64) 



which predicts the relative number of vector and pseudoscalar mesons. For the (cc), (be) and 
(bb) mesons the ratio (64) gives the following numbers: 0.46, 0.58, 0.52. The obtained results 
for the relativistic and bound state corrections to the different heavy quark fragmentation 
functions are shown in Figs. 2-5. Relative order contributions of relativistic and binding 
corrections are the biggest for the (be) mesons because of the little growth the parameter 
(60) and bound state energy W as compared to (cc) states. The decrease of these corrections 
in the bottomonium is explained by the increase of heavy quark mass (m c — > m&). All 
considered effects in the production of (be), (bb) and (cc) mesons are computed to the 
leading order in a s in the color singlet model. For S-states like J/^f, r] c , B c , B*, i] b and T 
the color-octet terms are suppressed relative to the color singlet terms by a factor of the 
second order over the relative velocity vq ^ . Our results should be useful for the comparison 
with more accurate (cc), (bb) and (be) meson production measurements in Z° decays or in 
pp collisions at the Tevatron. 
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